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Abstract. A combination ol measurements ol a s in e + e~ annihilation at LEP is presented. Distributions 
of various event-shape variables measured by ALEPH, DELPHI, L3 and OPAL at centre-of-mass energies 
from 41 GeV to 206 GeV are analysed in a common theoretical framework using 0(as)+NLLA predictions. 
The dominant theoretical uncertainties associated to missing higher orders are studied in detail. 



O 

OS 
O 

m 
O 

">< 

i 

Oh! 



X 



1 Introduction 

The most important parameter of QCD, the strong cou- 
pling constant a s , has been determined at LEP in different 
processes and with various observables over a wide range 
of energies. Here measurements using event-shape vari- 
ables provided by the four LEP experiments are combined. 
Six observables are analysed, these are thrust, heavy jet 
mass, wide and total jet broadening, C-parameter and 
the three-jet resolution parameter — In 2/3. The data have 
been collected at the LEPII-energies in the range from 133 
GeV to 206 GeV, at the Z boson peak and at lower ener- 
gies from 41 GeV to 85 GeV using radiative events. The 
combination of the measurements from different variables, 
energies and experiments takes into account correlations 
between the measurements. Theoretical uncertainties are 
re-evaluated for the individual input measurements and 
propagated to the combined values. This combination is 
performed by the LEP QCD working group PP . 



2 Input measurements 

The latest, partially still preliminary measurements of a s 
provided by ALEPH 0, DELPHI 0, L3 gj and OPAL [5] 
serve as input data for the combination. Altogether 194 
input values enter the combination precedure. The mea- 
surements are consitently extracted from the event-shape 
distribution using a theoretical description recommended 
by the LEP QCD working group pQ. Statistical and ex- 
perimental systematic uncertainties are evaluated by the 
experiments. In order to study hadronisation uncertainties 
the experiments have provided measurements using three 
generators PYTHIA, HERWIG and ARIADNE for hadro- 
nisation corrections. The dominant uncertainty is pertur- 
bative, and it is highly correlated between the input mea- 
surements. A special treatment is applied to calculate the 
theoretical uncertainty for each input measurement, as ex- 
plained in the next section. 



3 Theoretical predictions and uncertainties 
3.1 Predictions 

To second order in ct s , the distribution of a generic event- 
shape variable y (y = 1 — T, p, By/, Bt, 2/3, C) is given by: 
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tions are applied to the cumulative cross section 
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The fit function consists of two components, the fixed or- 
der term and an expression resumming leading and next- 
to-leading logarithms to all orders in a s . The recommended 
[H] prescription to merge the two calculations and to sub- 
tract off double counting terms is the Log(R) matching 
scheme 

lnR(y,a s ) = Lgi(a s L) + g 2 (a s L) - (G n L + Gi 2 L 2 )a s 
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with L = ]n(y /y) where y = 1 for y = l-T, p, y 3 , B T: B w 
and 2/0 = 6 for C. An improved version, the modified 
Log(R) scheme is used in the fit functions to determine 
a s , which vanishes at a given phase space limit y m ax- To 
fulfil this constraint L is replaced by 



L = - In 
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with the modification degree power p — 1 . A detailed col- 
lection of the formulae and numerical values for y max and 
the matching coefficients are given in 
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3.2 Uncertainties 



4 Combination procedure 



Theoretical uncertainties stemming from unknown higher 
orders in the perturbation series are inherently difficult 
to assess. The method adopted here includes a variation 
of the renormalisation scale x M from 0.5 to 2 and a new 
test which re-scales the logarithmic variable L — > In 

in the range 2/3 < it < 3/2, as suggested in 0. In addi- 
tion a different matching scheme and different implemen- 
tations of the kinematic modification are tested. Different 
sources of the perturbative uncertainty are combined with 
the uncertainty band method. This method calculates the 
uncertainty for a given variable and experimentally used 
fit range with a common value of a s (Mz), the latter being 
the result of a first combination iteration of all input mea- 
surements. The uncertainties of the distributions are cal- 
culated first, taking the prediction of the modified Log(R) 
matching scheme as reference theory and constructing an 
uncertainty band from the largest deviation with respect 
to the reference in each bin out of the theoretical variants 
mentioned above. In a second step the reference theory is 
calculated with variable a s and the smallest rcsp. largest 
value yielding a prediction inside the uncertainty band 
defines the theoretical uncertainty of a s . The method is 
illustrated in fig. ^ taking C-parameter as example. The 
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Fig. 1. Theoretical uncertainties for C-parameter. The grey 
area corresponds to the total perturbative uncertainty, dotted 
and dashed lines show individual contributions. The full lines 
represent reference predictions with values of a s shifted by its 
uncertainty. 



main two components to the uncertainty are the varia- 
tions of the scales x^ and xl , which contribute to a similar 
amount but in different regions of the distributions. 



The procedure to combine the 194 measurements begins 
with the construction of a 194 x 194 covariance matrix V, 
which relates the uncertainties for all pairs of input mea- 
surements. The covariance matrix is decomposed into four 
components corresponding to the four sources of errors: 
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The best combined value a s yielding the minimal error is 
a weighted average with weights obtained as follows 
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The uncertainty of the combination is given for each com- 
ponent as function of the weights, e.g. for the statistical 
uncertainty of tat = J2ij w i V^ at ' Wj. A specific treatment 
of the correlation of each uncertainty component concern- 
ing the off-diagonal matrix elements of V is applied. 

The statistical uncertainties are uncorrelated between 
different experiments and between different energies, but 
correlated for different observables using the same events. 
These correlations have been estimated numerically, using 
a large number of simulated datasets. 

The experimental systematic uncertainties are consid- 
ered to be uncorrelated between different experiments. Be- 
tween different energies and / or different observables of the 
same experiment, the "minimum overlap" assumption is 

Studies of the LEP QCD working group revealed that 
the hadronisation uncertainty is largely uncorrelated be- 
tween experiments because of different tunings and ver- 
sions of the generators. Therefore, the off-diagonal matrix 
elements of ]/ had - are set to zero and the hadronisation 
uncertainty is determined by repeating the combination 
procedure for each of the generators separately. The nom- 
inal results uses PYTHIA for corrections and as uncer- 
tainty the standard deviation of the three combinations is 
taken. For the combinations by energy the raw hadroni- 
sation uncertainty is fitted by a power law form a + b/Q 
in order to suppress statistical fluctuations. 

The theory uncertainties are expected to be highly cor- 
related between all pairs of measurements of the same ob- 
servable and to varying extents between measurements of 
different observables. The exact determination of the cor- 
relation, however, turned out to be difficult. Various at- 
tempts indicated correlation coefficients of the order of 
90%, but the result of the combination, in particular its 
theoretical uncertainty depends strongly on the assump- 
tions made. A more reliable approach is chosen, setting 
again the off-diagonal matrix elements of V th ' to zero and 
repeating the combination for two additional alternative 
input values of a° ±Z\a^ h ', where a* is obtained with the 
uncertainty band method. As perturbative uncertainty the 
difference between the nominal combined result and the 
alternative combinations is taken. 
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Table 1. Combinations of a s (Mz) for different energy ranges. 





Q <M Z 


Q = M Z 


Q> M z 


All 


a s {Mz) 


0.1208 


0.1200 


0.1201 


0.1201 


stat. error 


±0.0012 


±0.0002 


±0.0005 


±0.0003 


exp. error 


±0.0023 


±0.0008 


±0.0010 


±0.0009 


had. error 


±0.0032 


±0.0010 


±0.0007 


±0.0009 


th. error 


±0.0052 


±0.0048 


±0.0044 


±0.0046 




-0.0050 


-0.0048 


-0.0045 


-0.0047 


total error 


±0.0066 


±0.0050 


±0.0047 


±0.0048 



5 Results 

The combination is carried out for different subsets of vari- 
ables and energies, the full results are given in I]. The 
main result is summarised in Tabled where all variables 
and experiments are combined for different ranges in en- 
ergy. The best combination with the smallest total uncer- 
tainty is obtained with the LEPII data alone, combining 
all data yields a larger uncertainty. This is mainly due 
to the theoretical uncertainty, which is decreasing with 
increasing energy, scaling with a 3 s . The hadronisation un- 
certainty is deceasing as well, scaling with l/Q. 

The perturbative uncertainties are the key element 
limiting the precision of present elements, and their ex- 
act determination is cumbersome. As a cross check, re- 
sults from different variables are compared in fig. [2 As 
each variable receives possibly different higher order con- 
tribution, the spread of results in a s indicates their size. 
The RMS between results from different variables is 0.0016 
with a maximum spread of 0.0046, in good agreement with 
the perturbative uncertainty of ±0.0047 for the combina- 
tion of all data. 
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Fig. 2. Combinations of a s per variable. 



In order to study the energy evolution of a s (Q) the 
input measurements from all variables and experiments 



are combined per energy Q. The result is shown in fig. [21 
and compared to a the 3-loop evolution fit, which is found 
in excellent agreement with the data (x 2 /Ndo/=ll-7/13). 
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Fig. 3. Energy evolution of a 3 (Q). 



6 Conclusions 

A preliminary combination of measurements of a s from 
the LEP experiments using event-shape variables is pre- 
sented. The values of a s have consistently been deter- 
mined in a common theoretical framework and a new 
method is applied to assess perturbative uncertainties. 
The combined results are consistent with the expected 
energy evolution of a s . The preliminary result using all 
LEP data is a s (M z ) = 0.1201 ± 0.0048. 
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